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ABSTRACT 

We investigate physical properties of molecular clouds in disc galaxies with different 
morphology: a galaxy without prominent structure, a spiral barred galaxy and a galaxy 
with flocculent structure. Our 7V-body/hydrodynamical simulations take into account 
non-equilibrium H 2 and CO chemical kinetics, self-gravity, star formation and feedback 
processes. For the simulated galaxies the scaling relations of giant molecular clouds 
or so called Larson’s relations are studied for two types of a cloud definition (or 
extraction methods): the first one is based on total column density position-position 
(PP) datasets and the second one is indicated by the CO (1-0) line emission used 
position-position-velocity (PPV) data. We find that the cloud populations obtained 
by using both cloud extraction methods generally have similar physical parameters. 
Except that for the CO data the mass spectrum of clouds has a tail with low-massive 
objects M ~ 10^ — 10^ Mq . Varying column density threshold the power-law indices 
in the scaling relations are significantly changed. In contrast, the relations are invariant 
to CO brightness temperature threshold. Finally, we find that the mass spectra of 
clouds for the PPV data are almost insensitive to the galactic morphology, whereas 
the spectra for the PP data demonstrate significant variations. 

Key words: ISM: clouds - ISM: molecules - ISM: structure - stars: formation - 
galaxies: spiral - galaxies: structure. 


1 INTRODUCTION 

Molecular gas in galaxies is mostly concentrated in cold 
clouds with masses ~ 1 O‘*“®M 0 which are usually called gi¬ 
ant molecular clouds (GMCs). Their evolution is important 
for understanding transition of gaseous component into the 
stellar one. Indeed, galactic star formation generally occurs 
in the dense medium of GMCs. [Larson | |T98T| ) initially intro¬ 
duced three empirical scaling relations for the nearby molec¬ 
ular clouds in the Milky Way (MW). These relations reflects 
the general view on the GMCs properties and have the fol¬ 
lowing sense: 

• cloud size - line-of-sight velocity dispersion relation, 
(Tv oc , is the first one, it argues that the cloud structure 
is supported by the internal turbulence; 

• cloud virial mass - luminosity in CO lines relation. 
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Mvir OC Tqo) is the second one, it shows that GMCs are 
structures in the virial equilibrium; 

• luminosity in CO lines (sometimes cloud mass used) 
- size, Leo oc , is the third one, it claims that the mean 
cloud surface density Eo is likely to be constant if /Is « 2 . 


Despite a long way of the scaling relations investiga¬ 
tion, a complete theoretical explanation for the origin of the 
relations has not been offered yet. Based on the CO ob¬ 
servations of molecular clouds in the Galactic disc it has 
been found that GMCs have approximately constant sur¬ 
face density ~ 170 Mq pc“^ and the state of the clouds is 
really close to the virial equilibrium (Solomon et al.|[l987|. 


Roman-Duval et al. (20101 have found tight power-law cor¬ 


relation with index 2.36 ± 0.04 between radii and masses 
of the Galactic molecular clouds. The virial parameter of 
the derived clouds is mostly below 1 with the mean value 
0.46, so that clouds are strongly self-gravitating. Using ^^CO 
data Heyer et al. (20091 re-examined the scaling relations for 
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the Galactic clouds under assumption of constant CO-to-H 2 
conversion factor within a cloud. This leads to lower me¬ 
dian mass surface value, which is 42 M© pc“^. Note that 
the clouds found in this study are mostly unbound, that is 
in contradiction to the previous studies. Thus, the observa¬ 
tional data demonstrates signihcant scatter in the physical 
state of GMCs even in the Milky Way. 

For molecular clouds in both dwarf and giant disc 


galaxies Bolatto et al. (20081 have found scaling relations 
similar to that for the Milky Way clouds. They have con¬ 
cluded that GMGs identihed on its GO emission are the 
unique class of objects that exhibits a remarkably uniform 
set of properties from galaxy to galaxy. Meanwhile more 
recent comparison of GMGs in nearby galaxies by |Hughes| 
et al. ] poT^ let to hgure out that the GMGs properties 


(mass, radius, velocity dispersion) are not robust towards 
to the external conditions: clouds are smaller and fainter 
in less dense regions, i.e. inside low-mass galaxies and the 
outer regions of the Galaxy, compared to molecular struc¬ 
tures in denser environment, e.g. in the inner part of the 


Galaxy and other spirals like M 51 (Colombo et al. 20141 
and M 33 ( Engargiola et aL]|2003 Bigiel et al.|201(i I. 


Certainly, the scaling relations can reflect some univer¬ 
sality in both physical conditions inside clouds and inter¬ 
action of clouds with the ambient medium. Giant molec¬ 
ular clouds properties and evolution are governed by the 
interplay between self-gravity, magnetic field and feedback 
processes from stars born inside clouds. In many theoret¬ 
ical studies there has been attempted to understand how 
various feedback processes influence on the properties of 


GMCs (Shetty & Ostriker|2008 Taskerj 

201 1| Hopkins et al. 

2012 Braun et al.||2014l. For instance. 

Dobbs et al. (20111 


have traced the evolution of individual clouds in detail and 
found that cloud-cloud collisions and stellar feedback can 
regulate the internal velocity dispersion and lead to for¬ 
mation of unbound GMCs. Contrary to the previous study 


Tasker & Tan 120091 suggested that molecular clouds are 


gravitationally bound because of the low collisional rate of 
clouds relatively to its orbital time scale. Thus, the inter¬ 
nal turbulent energy can keep molecular clouds in the virial 


equilibrium. Several simulations of turbulence in GMCs (Re- 
|naud et al.||201^ [Kritsuk et al.|[2013[ ) have justihed that 
self-gravity plays an important role in the cloud structure, 
but doesn’t strongly affect on the ’velocity dispersion - size’ 
relation. 

Using high resolution simulations [Benincasa et ah] 


(20131 analyzed the physical properties of clouds whose 
number density is above 100 cm“®. They found that the 
slopes of the ’velocity dispersion - size’ and ’mass-size’ re¬ 
lations appear to be much steeper than the observational 
ones. On the other hand, [Tasker &; Tan| ( |2009| ) got a good 
agreement between the mass, radius, velocity dispersion of 
GMCs with those observed in the Galaxy. Such contrary 
conclusions are explained by not only differences in simu¬ 
lations, e.g. taking into account star formation and other 
processes, but also variety in samples of clouds caused by 
using different methods of cloud extraction. Moreover, |Fuji-| 


moto et al. (20141 found a significant effect of galactic envi¬ 


ronment on the cloud properties in the dynamical model of 
M 83. At first they established that the ’mass-size’ relation 
has bimodal distribution, and at second, GMCs tend to be 


less gravitationally bound in denser environment, i.e. spiral 
arms or bar, than in rarehed ones, e.g. inside disc. 

In numerical simulations a cloud is usually defined as 
an object whose gas density (column or volume) is higher 
than a given threshold. Such an object can consist of several 
dense molecular cloud lets surrounded by diffuse intercloud 
molecular and/or atomic gas. In addition there are several 
other methods for cloud definition based on dust extinction, 
molecular or/and atomic column density or CO intensity. 
For each method it is interesting to find the scaling rela¬ 
tions and compare it with the empirical ones established by 
Larson. That probably allows us to understand better what 
ISM structures are responsible for appearance of these rela¬ 
tions. 

The matching of the observed and simulated GMCs 
properties is not obvious because of different approaches 
used for cloud definition. In general, this problem has no 
unique solution both in observations, because in observa¬ 
tions the border of a cloud can depend on a chosen signal- 
to-noise limit. In numerical simulations there are two com¬ 
monly used methods for cloud extraction. The first one 
is based on total column density position-position (hence¬ 
forth PP) data sets and the second one is indicated by 
the CO line emission used position-position-velocity (PPV) 
data. The latter is utilized in CLUMPFIND (Williams et al. 
19941, CPROPS ( Rosolowsky fc Leroy|2006 l packages. 

In this paper we consider the physical properties ( 
namely, mass, radius, surface density, velocity dispersion, lu¬ 
minosity etc.) of clouds for two methods of cloud extraction 
based on PP and PPV datasets. In our simulations we study 
the scaling relations or so-called Larson’s laws for three MW- 
size galaxies with different morphology. The paper is orga¬ 
nized as follows. Section 1^ contains the description of our 
numerical model. Section [s] describes methods of cloud defi¬ 
nition. In Sectionj^we present the statistical analysis of the 
physical properties of molecular clouds. Section describes 
the scaling relations, the dependence of power-law indices of 
the relations on threshold value and mass spectra of GMGs 
for the simulated galaxies. In Section we summarize our 
key results. 


2 MODEL 

To simulate the galaxy evolution we use our code based 
on the unsplit TVD MUSGL (Total Variation Diminishing 
Multi Upstream Scheme for Conservation Laws) scheme for 
gas dynamics and the V-body method for stellar compo¬ 
nent dynamics. In gas dynamical approach we reach the 
second order in time and the third order in space using 
the minmod limiter. For the Riemann problem solution we 
adopt the HLLC (Harten-Lax-van Leer-Contact) method. 
More details about gas dynamic part of our code can be 


found in the paper Khoperskov et al. (20141. Stellar dynam¬ 


ics is calculated using the second order flip-flop integrator. 
For the total stellar-gaseous gravitational held calculation 
we use the TreeCode approach. 

For all models presented here we use a uniform grid 
with 4096 X 4096 x 512 cells for gas dynamics and set a 
computational domain 40 x 40 x 3 kpc with spatial resolution 
6 pc. The initial number of stellar particles is equal to 0.5 x 
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Table 1. Initial parameters adopted in the simulations. Here the following notations are assumed: is the mass of dark matter halo 

within 12 kpc sphere, is the halo scale length, crr(O) is the central radial velocity dispersion, cTz/c^’r is the ratio of the vertical velocity 
dispersion to the radial one, Sgo is the central surface density of gaseous disc, hg is the radial scale length of the gaseous disc, A^tot 
the number of clouds extracted with the tiireshold (CDN is a cloud definition when the total column density within the cloud 

exceeds the AT^ot threshold and the cloud is delineated by corresponding level of the total column density), A^cf is the number of clouds 
extracted using CLUMPFIND (henceforth we use CF abbreviation for shortness). The following parameters are the same for considered 
models: S*o = 835 Mq pc“^ is the central stellar surface density, A* = 3 kpc is the radial scale length of stellar disc, and ha^ = 2h* is 
the radial velocity scale length, Sgo = 10 Mq pc~^ is the central gas surface density. 


Model (Morphology) Halo Bulge Stellar disc Cloud definition 

Mh tth Mb 6b CTz/^r AfcDN A^CF 

10^^ Mq kpc 10^^ Mq kpc km 


B (No structure) 

8.8 

3.857 

- 

- 

75 

0.5 

1095 

1150 

F (Milky Way like) 

GO 

00 

1.1 

0.7 

0.153 

100 

0.7 

1065 

1203 

H (Flocculent) 

8.25 

1.1 

- 

- 

50 

0.45 

1012 

1111 


10®, during the simulation it reaches 2 x 10® depending on 
star formation activity. 


2.1 Chemical kinetics and gas thermodynamics 

Usually the emission in CO lines is the major source of the 


information about the GMCs (Dame et al. 2001 Bolatto 
|et al.|2008||Leroy et al.|2009[ ), and the intensity in CO lines 
is used to restore the mass of molecular hydrogen through 
Xco factor ( Dickman||1975 Bolatto et al.|[2013 |. Then, we 
are interested in a reasonable CO chemical network that on 
one hand gives fine CO molecule evolution and on the other 
requires adequate computational resources. Rather detailed 
networks include more than 20 chemical species involved 
in several hundreds of reactions (e.g. Omukai||2000 Glover 


et al.|[20l0 l, which is computationally unacceptable for our 


purposes. Fortunately, Glover & Clark (20121 found that 


the reduced network proposed by Nelson & Langer (19991 


gives adequate results in comparison to the detailed chemical 
model, which consists of 218 reactions amongst 32 species 
I Glover et al.|20l0 1 . So that here we exploit the model based 
on the network proposed by Nelson & Langer (19991. 


Based on our simple model for H2 chemical kinet¬ 
ics I Khoperskov et al.|2013 1 we expand the Nelson & Langer 


119991 network by several reactions needed for hydrogen ion¬ 


ization and recombination. For H2 and CO photodissocia¬ 
tion we use the approach described by [Draine fc Bertoldil 
119961. The CO photodissociation cross section is taken from 


Visser et al. (20091. In our radiation transfer calculation 


described in Section 2.3 below we get ionizing flux at the 
surface of a computational cell. To calculate self-shielding 
factors for CO and H2 photodissociation rates and dust ab¬ 
sorption factor for a given cell we use local number densities 
of gas and molecules, e.g. = nu^L, where nH2 is H2 
number density in a given cell and L is its physical size. 
The chemical network equations is solved by the CVODE 
package ( Hindmarsh et al.|2005 |. 

We assume that a gas has solar metallicity with the 
abundances given in|Asplund et al.| (|2005|): [C/H] = 2.45 x 
lO”"^, [0/H] = 4.57 X 10"'*, [Si/H] = 3.24 x 10“^ Dust de¬ 
pletion factors are equal to 0.72, 046 and 0.2 for C, O and 
Si, correspondingly. We suppose that silicon is singly ionized 
and oxygen stays neutral. 

For cooling and heating processes we extend our previ¬ 


ous model ( Khoperskov et al.|2013 l by CO and OH cooling 
rates (Hollenbach & McKee 19791 and Cl fine structure cool¬ 
ing rate ( Hollenbach fc McKee|19^ I. The other cooling and 
heating rates are presented in detail in Table 2 (Appendix 
B in Khoperskov et al.|2013 1. Here we simply provide a list 
of it: cooling due to recombination and collisional excita¬ 
tion and free-free emission of hydrogen ( Cen||1992 1, molec¬ 
ular hydrogen cooling (Galli & Falla 19981, cooling in the 


fine structure and metastable transitions of carbon, oxygen 
and silicon ( Hollenbach fc McKee|[r989 ), energy transfer in 
collisions with the dust particles (Wolfire et al. 20031 and 
recombination cooling on the dust ( Bakes fc Tielens||1994 1, 
photoelectric heating on the dust particles (Bakes & Tielens 
T9MI [Wolfire et al.||2003| , heat ng due to H2 formation on 
the dust particles, and the H2 photodissociation (Hollenbach 
fc McK^|1979 1 and the ionization heating by cosmic rays 
( Goldsmith fc Langer|[T978 1. In our simulations we achieve 
gas temperature value as low as 10 K and number density 
as high as 5 X 10® cm“®. 


2.2 Star formation and feedback 

In the star formation recipe adopted in our model mass, 
energy and momentum from the gaseous cells, where a star 
formation criterion is satisfied, are transited directly to new¬ 
born stellar particles. A star particle is formed in a grid cell, 
if the following criteria are fulfilled: (i) the gas density in the 
cell should be higher than 100 cm“® (such a high value pre¬ 
vents the formation of huge number of stellar particles), (ii) 
the total mass of gas in surrounding cells exceeds the Jeans 
mass Mceii > Mj (this help us to avoid the star formation in 
hot and warm medium, where some feedback processes oc¬ 
cur), and we adopt the local star formation efficiency 0.01. 
In star-forming cells the number density and temperature 
reach n > 200 — 500 cm and T ~ 50 K. 

Feedback model includes several sources of thermal en¬ 
ergy, namely stellar radiation, stellar winds from massive 
stars and SN explosions. The amount of injected energy con¬ 
nected to these processes is calculated for each stellar par¬ 
ticle using the stellar evolution code STARBURST99 ( jLei-j 
therer et al. 19991. We model supernova feedback only as 


thermal energy injection into a gas. We take into account 
mass loss by stellar particles due to SN explosions and stel¬ 
lar winds from both massive and low-mass stars. 
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Figure 1. Initial conditions for simulated stellar-gaseous discs. Circular velocity is shown in the left panel. Stability parameter Qy 
adopted for finite thickness two component disc and using the approximation by [Romeo fc Wiegert| ( [2011^ is shown in the right panel. 
Physical parameters of galaxy models are presented in Table [T] 


2.3 Radiation transfer 


To account molecule photodestruction we should know spa¬ 
tial structure of UV background in the galactic disc. Recent 
observations provide some evidences for significant radial 
and azimuthal variations of UV flux in the nearby galax¬ 
ies (Gil de Paz et al. 2007). No doubt that such variations 
are stipulated by local star formation. So that we need to 
include radiation feedback from stellar particles in our cal¬ 
culations. 

Through our simulations the UV emission of each stellar 
particle is computed with the stellar evolution code STAR- 
BURST99 ( Leitherer et al.||1999 l assuming solar metallic- 
ity of stellar population. So that for each particle we know 
its luminosity evolution. After that we separate particles in 
two groups: young stellar particles (the age is smaller than 
20 Myr) and the other ones. For definiteness we assume a 
uniform background field ten times lower than that in the 
Solar neighbourhood, F;, = 0.1 Habing.Thus the UV back¬ 
ground in a hydrodynamical cell with coordinates ro 
can be written as 


F'^'^(ro) = Fb -f ^ (ro) + ^ r,), (1) 

i 3 

X]iF°*‘*(ro) is deposit from old stellar population (age > 
20 Myr), which plays a role only in a cell where the stel¬ 
lar particle locates (ro). The last term is UV flux from 
young stellar population - the brightest stars. Their deposit 
is the most important in photodestruction of molecules in 
surrounding medium. 

Due to number of young stars is small at each time step, 
we can use the ray-tracing approach for each stellar particle. 
For j-th ’’young particle” we estimate the radius of spherical 
shell (similar to the Stroemgren sphere), where the UV field 
value decreases down to 0.1 Habing: 

F," = 0.iyF*/(4^), (2) 

where L* is luminosity of y’-th stellar particle in Habing units 
and 5 = \/ (Sx)^ + (Sy)^ + (Sz)^ is effective cell size. For 
each shell we calculate the UV flux assuming the optical 
depth r = 2V/(10^^cm“^), where N is the total column 
density of gas in cm“^. So that we can get the distribution 


of the UV intensity in the entire galactic disc according to 

Eq.[T] 


2.4 Model of galaxies 


We start our simulations from the self-consistent radial and 
vertical equilibrium state of stellar-gaseous discs in the fixed 
gravitational potential of dark matter halo. We assume both 
stellar and gaseous discs have exponential form, but with dif¬ 
ferent spatial scale lengths. Circular velocity of the gaseous 
disc embedded into the gravitational potential can be found 
as: 

( Q^halo d^'bulge d^'diac d^'gaa \ , . 

r \ dr dr dr dr ) ' 


where 4Ihaio is gravitational potential of dark matter halo, 
the halo is assumed to be steady isothermal sphere, Tbuige 
is potential of bulge, Tdisc is potential of stellar disc and 
Tgaa is potential of gas. The parameters of the gravitational 
potential can be found in Table 0 Fig- a presents the ra¬ 
dial dependence of circular velocity for the galactic models 
considered here. 

For stellar particle kinematics the asymmetric drift is 
taken in the form of Jeans approximation: 


- a" 


/ r 1 d(E.a?) \ 

V 0-2 E* o-p dr J ’ 


(4) 


where Ur, o^p are radial and azimuthal velocity dispersions, 
respectively, E* is stellar surface density distribution. The 
parameters of the potential and matter distributions can be 
found in Tabled To compute an initial distribution of stars 
we solve Eg. [4pus ing the iterative procedure described in 
Khoperskov etah] ( |2003 K 

Despite the parameters of the models of galaxies pre¬ 
sented in Table are close to each other, the various sta¬ 
bility conditions allow us to follow galaxies with the differ¬ 
ent morphology. Initial stability criteria for two-component 
models (stellar-gaseous) are shown in the right panel of 
Fig.0 We compute three models of the stellar disc equi¬ 
libria: gravitationally over-stable disc (Model F, without 
prominent structure), highly unstable (Model H, flocculent 
spirals morphology) and intermediate state disc (Model B, 
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Figure 2. The projected maps (inclination 30°) of stellar surface density (top left group of panels), stellar UV radiation field (top right 
group of panels), total gas column density (left group of panels in the middle row), radial (line-of-sight) velocity component (right group 
of panels in the middle row) and CO integrated intensity (bottom group of panels) at t = 500 Myr for the following models of galaxies: 
no spiral structure or model B (left map in all groups of panels), Milky Way-like or model F (central map in all groups of panels) and 
flocculent galaxy or model H (right map in all groups of panels). Initial parameters of the galactic models can be found in Table ^ 


MW-like morphology). The initial stability parameter for 
two-component disc model accounting the finite disc thick¬ 
ness effect is adopted in the form by |Romeo & Wiegert| 

i Mit (Fig.g. 

Figure shows the maps of the stellar surface density, 
stellar UV radiation field, total gas column density and CO 
integrated intensity at t = 500 Myr for the following mod¬ 


els of galaxies inclined by i = 30°: a galaxy withoutspiral 
structure (model B), a Milky Way type galaxy (model F) 
and a flocculent galaxy (model H). The initial parameters 
of the models are given in Table Note that the spatially 
averaged UV radiation field in all three models of galaxies is 
significantly greater than a value of 0.1 Habing (see upper 
middle row in the Fig.j^, so that our choice of the uniform 
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background is reasonable. We have adopted inclination angle 
i = 30° as a value which is enough to get signihcant line-of- 
sight velocity scatter while the structures in the gaseous 
disk are still rather well spatially distinguishable. 


3 CLOUDS DEFINITION 


Prior to the calculation of physical parameters of clouds, 
we should define what the cloud is. In the most obvious 
approach for cloud definition (CD) a cloud is an isolated 
gaseous clump with gas density (column or volume) higher 
than a given level. This is the simplest criterion, but it 
doesn’t reflect the chemical composition of a clump and we 
cannot say anything about molecular content of such cloud. 
Moreover the methods based on total gaseous column den¬ 
sity are not relevant to observable values, because using such 
methods some material, which is not associated with cloud 
itself and laid along the line of sight, can be regarded as 
a part of a cloud. This reveals in physical parameters of a 
cloud. So that we need a criterion based on the distribution 
of molecules in the ISM. Such criterion connects both chem¬ 
ical and extinction properties of a cloud and allows us to 
separate two phases of the cold interstellar medium: atomic 
and molecular gas. 

Usually molecular clouds are stud ied by its emission in 
molecular lines (e.g., for (se e e.g. Solomon et al. 19871, 
(see, e.g. 


Heyer et al.|2009| and more recently in OH 


by Allen et al. (20151). Since our model includes the H 2 and 


CO molecule kinetics, we can use CD criteria based on both 
total gas column density and intensity in CO lines. This 
allows us to check the range of applicability for each CD 
criterion. So that the properties of a particular cloud are 
expected to depend significantly on the extraction criterion. 
It is also not clear how a choice of criterion influences on 
the statistical properties of the whole ensemble of molecular 
clouds. 

Below we consider two approaches related to the prop¬ 
erties of a cloud measured in observations. We dehne a cloud 
as a region inside that the total (molecular and atomic) 
hydrogen column density is higher than a given threshold 
A^tot (henceforth we use the abbreviation CDN for the 
method and corresponding indices). According to the CDN 
criterion we firstly find all local maxima (peaks) of the gas 
column density in the plane of the galactic disc. After that 
around each local maximum we find cells with value higher 
than a given threshold. In some cases these regions merge 
into larger ones and form a cloud with several local maxima. 
Usually such coalescences take place in dense galactic struc¬ 
tures, e.g. in spiral arms, bar or regions nearby the galactic 
center. Thus, our approach for finding clouds is a combi¬ 


nation of the ’contour method’ (Fujimoto et al. 20141 and 
’peaks method’ | Tasker fc Tan|[2009 1. 

One of the widely used method for extracting struc¬ 
tures from PPV data cubes is CLUMPFIND ([Williams et al.| 


1994). This method is based on contouring data array at 


many different levels starting from the peak value and mov¬ 
ing down to specific threshold. In the present work the CU¬ 
PID implementation of the CLUMPFIND algorithm is used 


for CO intensity method of the clouds extraction (Berry 
jet al.|2007| ). Henceforth we use the abbreviation CF for this 
method and corresponding indices of variables. 


We calculate CO brightness temperature in the form of 
the PPV data cubes using the method described in [Feld-j 


mann et al. (2012). In the calculations the spectral velocity 
resolution equals to 0.5 km s“^ , which is potentially enough 
to resolve structure of massive clouds.This spectral resolu¬ 
tion is comparable to that one reached in the recent inter¬ 
ferometric observations (see e.g., Roman-Duval et al.|2009 |. 
Note that we discuss the dependence of power-law indices of 
the relations on velocity resolution value in the Section 5.4. 

Certainly, using two above-mentioned criteria we get 
two different population of clouds. Number and total mass 
of clouds are also different and depend on the value of col¬ 
umn density and brightness temperature thresholds. In our 
analysis we take = 1.9 x cm“^ and = 3 K 
as fiducial threshold values, which provide us comparable 
numbers of clouds (around 1000) and similar total gaseous 
masses locked in clouds (Mt « 2 — 4 x 10® Mq ). Note that 
these values of Mt are close to the total mass of molecular 
clouds in the Milky Way ( [Williams fc McKee||19f)7| ). The 
numbers of clouds and total gaseous masses for the galactic 
models considered here are given in Table 

For instance, a small region with the spatial distribu¬ 
tion of extracted clouds in the MW type galaxy (model F) 
is shown in the Fig. We mark the extracted clouds as 
coloured areas in contrast to the grey scale background of 
gaseous column density map. It is clearly seen that spatial 
distributions and numbers of clouds are remarkably dis¬ 
tinct for the considered criteria. Prior to any quantitative 
analysis of the physical parameters we should notice two 
issues. On one hand, non-interacting clouds are appeared 
to have similar shape for both extraction methods. How¬ 
ever, in more dense environment clouds extracted by differ¬ 
ent methods look very unalikely. We suppose that this can 
be a result of dynamical effects related to cloud collisions 
and/or stellar feedback effects. On the other hand, it seems 
that very large clouds (or agglomerations) extracted by us¬ 
ing the CDN method have internal structure, which we can 
hardly resolve because of our spatial resolution is still not 
enough high. Thus, for the CDN criterion large clouds and 
cloud chains (at least in the dense environment) can be ex¬ 
tracted, while using the CF method such large structures 
are splitted into individual lumps with internal motions and 
other specific inhomogeneities. 


4 GMCS PHYSICAL PARAMETERS 

Cloud formation were studied numerically in detail by 
[Dobbs et al.[ ( (2006|[2008| ); [Dobbs] ( [2008[ ). We mention that in 
our simulations clouds are result of self-gravity, thermal in¬ 
stability, cloud collisions and other processes occurred in the 
galactic disc. Here we briefly describe physical parameters 
of the cloud samples obtained in our analysis. 

On the one side, spiral arms stimulate GMCs formation 
due to gas falling into the gravitational potential well of the 
arms. Gravitational potential of spiral structure induces 
collisions of clouds that in turn stimulates star formation. 
From the opposite side, supernovae explosions in star form¬ 
ing spiral arms can destroy clouds. So one can conclude that 
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Figure 3. A region with the spatial distribution of extracted 
clouds in the model F. The extracted clouds are marked by 
color: top panel shows the distribution for CF (or CLUMPFIND) 
method, bottom panel demonstrates that for CDN approach. The 
background greyscale map is the gas surface density. 


of clouds in the inter-arm regions requires longer lifetime 
i Scoville et al.]1979 |Koda fc et al.|2()09 1. So that a question 
about lifetime of molecular clouds is still under debates (see 


e.g. Dobbs &l Pringle|2013 Zasov fc Kasparova|2014 1. 

Number of clouds extracted by using both criteria for 
the cloud definition described above depends on a choice 


1.9 X 10^^ cm~^ and = 3 K we extract ~ 1000 isolated 
clouds in our simulated galaxies. These clouds have the fol¬ 
lowing physical parameters: masses are « 10"^ — 10^ Mq , 
sizes vary within the range 3 — 100 pc, one dimensional ve¬ 


locity dispersion is ranged 0.1 — 10 km s“ , mean surface 
densities are ~ 60 — 300 Mq pc“^, and luminosities in CO 
lines are 10^ — 10^ K km s”'^ pc^ . These parameters de¬ 
pend slightly on the galactic morphology. Figure shows 
the distributions of these physical parameters for the cloud 
population in the Milky Way type galaxy (Model F) for both 
techniques of the cloud definition. 

For the CF method the physical paramters of clouds, 
i.e. masses, ID line-of-sight velocity dispersion (LOSVD), 
total luminosity and cloud sizes, are calculated using the 


prescriptions from Williams et al. (19941. 


For the CDN approach one dimensional velocity disper¬ 
sion of a cloud is calculated according to 




(5) 


where Uc is cloud center mass velocity vector, u is cloud ve¬ 
locity vector. Such approach is widely used to extract clouds 
in numerical simulations (see e.g. Benincasa et al.|2013 Fu- 
jimoto et al.||2014 and references therein). A cloud size is 

[a 

calculated according to R — \ —, where A is cloud surface 

V TT 

in pc^. 

The ratio between the kinetic energy and gravitational 
energy is commonly used to specify deviation from the virial 
state of a cloud under assumption of the constant density 
distribution ([Bertoldi fc McKee[|1992|: 




GMiu 


Mlum 


( 6 ) 


where a is cloud virial parameter, Rd is cloud size in parsecs, 
Mlum is mass of cloud in Solar units adopted from its CO 
luminosity using Xco conversion factor ( Dickman||1978 (: 

4.4LcoAco 


Mlum = 


2 X 1020 


(7) 


It is easy to see that the luminosity mass Mium and the virial 
mass Mvir of a cloud are generally not equal to each other: 

Mvir ~ aMlum • (8) 

Middle panels in Fig. present the parameter a for cloud 
population in the Milky Way type galaxy. It is clearly seen 
that the majority of molecular clouds tend to be in the 
virial equilibrium. Note that the distribution of such quasi- 
virialized objects is close to the uniform distribution as a 
function of cloud mass. The physical parameters obtained 
for our models of galaxies are in agreement with other re¬ 


cent numerical simulations (Tasker & Tan 2009 Tasker 2011 
Khoperskov et al.|2013|. Such result is likely to be a reflec¬ 


tion of the turbulent energy distribution in the entire galac¬ 
tic disk (detailed study see in Kraljic et al.||20i4 (. 

The distributions of cloud masses obtained by using the 
CF and CDN methods are slightly different: in the former 
we extract smaller and less massive clouds than the ones 
extracted in the latter |^,b). Moreover, the mass range for 
the CF sample of clouds is rather wide. The reason of this 
is clearly seen in the Fig.[^ large structure extracted by the 
CDN method are divided into several smaller clouds when 
the CF technique is used (see Fig. |^. We discuss several 
dynamical and methodological effects related to this issue 
in the further paragraphs. 

The ID velocity dispersion of clouds defined by Eq. [^is 
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unlikely to provide a good description for observed line-of- 
sight velocity dispersion, making it higher at least for extra- 
galactic GMCs. Moreover, the regular quasi-circular motion 
of giant clouds around galactic center leads to overestima¬ 
tion of velocity field within a cloud about 1-2 km s“^ . But 
this effect is significantly smaller than the ID LOSVD value. 
Spectral resolution 0.5 km s“'^ allows us to distinguish inter¬ 
nal structure of large clouds and measure the ID LOSVD 
within 1-2 km s“^ accuracy (see Fig. [^). Using the CF 
method we extract more homogeneous cloud sample which 
has smoother (without many local peaks) distribution. In 
any case both methods provide more or less similar shape 
of the velocity dispersion distribution functions, which are 


close to the observable ones (see e.g. Roman-Duval et al. 


2010). It seems that the CF method splits large clouds into 


smaller ones due to the complex velocity structure, which 
mainly takes place in colliding and tangent gaseous flows or, 
in general, turbulent regions. Note a remarkable difference 
between clouds extracted by the CF and CDN methods: 
large clouds found by the CF method are isolate lumps lo¬ 
cated in calm environment, whereas large clouds extracted 
by using the CDN criterion mostly represent dynamically 
interacting structures. 


5 SCALING RELATIONS ANALYSIS 

In this section we discuss the scaling relations for CMCs 
extracted according to two criteria of the cloud definition 
for the galaxy models described above. In the Table we 
collect all indices and normalizations for the scaling relations 
obtained in the analysis of the models of galaxies considered 
here. It seems that there is no strong variation of the CMCs 
scaling relations obtained by using the CDN method for 
galaxies with different morphology. This is considered in 
Sects. |5.H |5.2[ |5.3| in detail. A role of the environment in 
the galactic disk on the CMCs parameters is discussed in 
the Sect. 15.61 

The statistical relations for the three models of galaxies 
considered here are presented in Figs. and described 

in the corresponding subsection. The top row of the panels 
in each figure shows the relations obtained by using the CF 
method, and the bottom one shows the relations based on 
the CDN criterion. 


largest dimension of the cloud. So the use of the CDN cri¬ 
terion can result in incorrect estimation of cloud sizes and 
overestimation of their column density and velocity disper¬ 
sion (Fig. 1^. 

From the bottom panels in Fig.|^one can conclude that 
for the adopted here column density threshold the clouds 
extracted by using the CDN criterion hold gas with higher 
velocity located at their periphery. Such intercloud (diffuse) 
gas can contain significant molecular fraction. Note that in 
some recent observations extended structures with signifi¬ 
cant molecular fraction are found around molecular clouds 
i Caldu-Primo et al.||2015 |. The velocity dispersion of these 
structures is higher than that in the clouds. This can be 
considered as an evidence that molecular gas can exist in 
two phases: clumpy phase, which is organized in molecular 
clouds, and diffuse one, which is located in extended struc¬ 
tures around clouds. 

CO molecules are efficiently formed only in the dense 
shielded environment and are destroyed due to heating and 
photodissociation by the stellar feedback. Using the CDN 
criterion total hydrogen column density has a deposit from 
not only central dense molecular regions, but also peripheral 
parts of a cloud that can mainly contain atomic hydrogen 
and even some intercloud starforming regions, where young 
stars have already existed. If we extract clumps brightly 
emitting in CO lines, then the low density HI gas at the 
periphery is excluded from the consideration. 

One can see that the indices in the power-law relation 
Cv — Rci for the models of galaxies with more pronounced 
structures significantly differ from that for the observational 


fits ( Larson|1981 Solomon et al.p987 Bolatto et al.|2008 |. 
This deviation takes place for both threshold criteria, but it 
is smaller for the CF method. This can be explained by the 
fact that using the CDN criterion we can consider gaseous 
structures which are not really associated with clouds. So 
that gaseous flows at the outskirts of cloud are added to the 
internal turbulence motions of this cloud, hence, the nu¬ 
merical ratio between velocity dispersion and size of cloud 
becomes higher. We mentioned above that the CF approach 
is a sharper ’filter’ for molecular clouds than the CDN one, 
and our cloud samples based on the CO datacubes data 
demonstrate the statistical relations closer to the observed 
ones. Thus, in our simulations the first Larson’s scaling re¬ 
lation is better reproduced for PPV data. 


5.1 Velocity dispersion - size relation 

Fig.|^ shows that the clouds extracted according to the CDN 
criterion have higher velocity dispersion with the mean value 
~8 —20kms“'^ compared to that obtained for the CF one, 
in this case the mean value of velocity dispersion decreases 
to ~ 1 — 5 km s“^ . One can note that the observational fits 
galaxy and others ( [Solomon et al.|[l987[ 
I are in better agreement with the rela¬ 
tions obtained for the CF sample of clouds. The relations for 
the clouds extracted by the CDN method show significant 
deviation from the observational fits. 

Clouds have extremely complicated shape and consist 
of crossed and elongated structures (see Fig.|^. So that high 
total hydrogen column density at the periphery of clouds can 
be a geometrical effect when the line-of-sight goes along the 


tor the Milky Way 


Bolatto et aI.|20TO 


5.2 Virial mass - luminosity relation 

The ’virial mass - luminosity relation’ reflects a suggestion 
that CMCs state is close to the virial equilibrium. Figui 
shows the correlation between virial mass Myir (see Eq 
and total luminosity of the extracted molecular clouds for 
three models of galaxies. One can see significant scatter 
of the physical parameters for the simulated clouds around 
the observational fits. Similar to the ’velocity dispersion 
- size’ relation one can see here also that the scatter for 
the CDN criterion is larger than that for the CF one, es¬ 
pecially this is remarkable for low surface luminosity clouds 
(Leo ~ 10^ K km pc^ ). In general, for the same lumi¬ 
nosity value the virial mass of clouds obtained for the CDN 
criterion is systematically greater than that for the CF one. 
That can be explained by that high virial masses have large 
clouds formed in collisions of smaller ones, such massive 
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Figure 4. The physical parameters of GMCs obtained using CDN criterion (blue) and CLUMPFIND algorithm (green) for the model 
F galaxy (middle panel of Fig.[^: size (a), mass (b), virial parameter (d), velocity dispersion (e) and luminosity (f). The dependence of 
virial parameter on mass of a cloud is shown in panel (c). The red dotted and black solid lines depicted in panels (c) and (d) correspond 
to a = 1 and a = 2, respectively. 


clouds are mainly associated with spiral arms and/or bar. 
During collision of clouds, molecules can be destroyed, but 
the shock waves cannot ionize gas (or such gas recombines 
rapidly). So that significant part (on mass) of such clouds 
is locked in the warm atomic hydrogen phase. Then, using 
the CDN criterion we obtain clouds with high total hydro¬ 
gen column density, where the deposit of atomic hydrogen 
to the column density is dominant or very significant. 

The CO brightest clouds are really molecular ones and 
they are probably to belong to older cloud population: a gas 


inside them has to become molecular (Glover et al. 20101. 
Whereas the darkest (massive) clouds in CO line are be¬ 
lieved to be either young population of massive clouds, in 
which atomic hydrogen hasn’t yet transformed into molec¬ 
ular, or maybe pseudo-virialized structures, which consist 
of a group of small molecular clouds ’bounded’ by atomic 
intercloud and/or more diffuse molecular medium (see|^. 
Such structures can appear in dense environment, e.g. spiral 
arms, bar, central parts of a galaxy, where a chosen thresh¬ 
old is low enough to extract separate clouds and leads to 
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Table 2. The parameters of the scaling relations for GMCs in the simulated galaxies. 


Model (Morphology) - CD 

(Tv = Al 
Al 

km s“^ pc“l^i 

^cl 

dl 

Mvi, = Aaiff 
A 2 

M 0 (K km s“^ pc^ )-/32 

(92 

Lci = AiR^J 
As 

K km pc^ pc~^^ 

ds 

B (No structure) - CDN 

0.94 ±0.25 

0.66 ±0.17 

795 ± 151 

0.64 ± 0.31 

200 ± 45 

1.72 ±0.13 

B (No structure) - CF 

0.97 ± 0.12 

0.49 ±0.1 

15.8 ±2.2 

1.56 ±0.11 

1000 ± 38 

1.02 ± 0.2 

F (Milky Way like) - CDN 

0.35 ±0.21 

1.11 ± 0.18 

16 ±8.1 

1.04 ± 0.28 

202 ± 29 

1.73 ±0.11 

F (Milky Way like) - CF 

0.57 ± 0.18 

0.69 ± 0.12 

13.1 ± 3.6 

1.47 ± 0.21 

630 ± 44 

1.28 ±0.19 

H (Flocculent) - CDN 

0.87 ± 0.24 

0.76 ±0.16 

1584 ± 212 

0.62 ±0.29 

156 ±51 

1.68 ±0.3 

H (Flocculent) - CF 

0.84 ± 0.9 

0.54 ± 0.09 

16.2 ± 3.1 

1.56 ±0.1 

1000 ± 35 

1 ± 0.21 
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Figure 5. The Velocity dispersion - cloud size’ relation obtained for the CF algorithm (top row of panels) and the CDN criterion 
(bottom row of panels). Left column of panels corresponds to the model of galaxy without spiral structure (model B), middle column 
presents the relation for the MW-like galaxy model (model F) and the right column shows the relation for galaxy with flocculent structure 
(model H, see Fig.[^. The solid red line is a power-law fit for the simulated data (the corresponding formula is shown in the right bottom 
corner). The dashed blue line corresponds to the fit for the data taken from |Bolatto et al.| | |20d^ . The dotted black line shows the fit to 
the data on the Milky Way clouds got by [Solomon et al.| ||1987|l. 


merger small clouds into larger structure. The use of the CF 
criterion provides more reasonable cloud sample and doesn’t 
lead to extracting such large gaseous structures. So that the 
scatter for the sample of clouds obtained for the CF criterion 
is much smaller that for the CDN one. 

For the CDN criterion the slope of the fit for the sim¬ 
ulated sample of clouds is flatter than that for the observa¬ 
tional data (see Fig.|^. Obviously, it comes from the excess 
of massive clouds with low surface luminosity. Using the CF 
approach the picture for all three models of galaxies shows 
opposite behaviour. The slope becomes steeper than that 
obtained in the observations. Above one can see that the 
CF method usually leads to splitting large structures into 
smaller ones due to systematic LOS velocity variations for 
a given structure (see Fig. [^, that reveals in both size and 
mass distributions (see Fig. W. Thus, we have a relatively 
large subsample of small clouds, which cannot be resolved 


(in space and/or in the LOS velocity coordinate) in obser¬ 
vational data. 

One can suppose that if a large gaseous agglomeration 
in the vicinity of spiral arms or galactic centre is split- 
ted into several isolated clouds, then the number of massive 
and bright clouds becomes lower, whereas smaller clouds are 
more numerous. Thus, the slope of the fit can become flat¬ 
ter. So that the increase of threshold value is likely to result 
in better match with observations. 


5.3 Luminosity - size relation 

Originally Larson (1981) found the ’mass-size’ relation for 
the Galactic molecular clouds: Md oc i?d- This relation can 
be interpreted as molecular clouds have the same (constant) 
surface density. Here we use another form of this relation, 
namely ’luminosity - size’, because it includes at least one 
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Figure 6. The ’cloud virial mass - luminosity’ relation. The other notations are the same as in the Fig.[^ 



Figure 7. The ’luminosity - cloud size’ relation. The dotted lines correspond to the cloud surface density according to Eq. for the 
constant conversion factor Xqq = 2 x 10^*^ cm“^ (K km s“^)“^ . The other notations are the same as in Fig. [s] 


observable value. Note that mass of a cloud can be easily 
found from luminosity using conversion factor Xco accord¬ 
ing to Eq. such re-calculation doesn’t affect on the slope 
of the scaling relation in case of the constant conversion fac¬ 
tor. 

Figurej^shows the relation for three models of galaxies. 
It is clearly seen that for all models the surface density of 
clouds is locked within interval ~ 10 — 1000 M© pc“^. This 
surface density range is rather universal within M = 10® — 


10^ Mq for all galaxy models. For the CDN criterion one can 
note substantial scatter of the cloud parameters below the 
critical value of the surface density « 10 ^ Mq pc“^ (see bot¬ 
tom row in Fig.j^. This is just a reflection of existence of the 
dimmer parts of the clouds. It seems that the strong limit 
on the maximum value of the surface density can be inter¬ 
preted as a result of ongoing star formation, which prevents 
the formation of more dense clouds. Molecules in such clouds 
are destroyed immediately due to photodissociation by UV 
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Table 3. The power-law indices of the scaling relations found in 
the Model F for several values of CO spectral line resolution 5v. 
The observational fits are given in the top frame of Table. 


Model / Observations 6v ^2 ^3 


km s 


Larson 


a98i|i 


Bolatto et al. 

to 

o 

o 

00 

yolomon et al. 

■ITM7 


Homan-Uuval et al. 20101 


0.38 

0.5 

0.6 


0.81 

1 


2.55 

2.54 

2.36 


Model F 

0.5 

0.69 

1.47 

1.28 

Model F 

1 

0.65 

1.34 

1.74 

Model F 

5 

0.6 

1.19 

2.41 


radiation from newborn stars. However, such picture cannot 
be supported by the analysis of the clouds extracted by CF 
method (see top row in Fig.[^. It is possible that the bright¬ 
ness of large clouds becomes lower than that expected due 
to the shielding effects. Note that the optical depth effects 
become important when the value of the column density ex- 


ceeds ~ 2 x 10^'^(r/10^) ^ cm ^ (e.g., Hollenbach & McKee 


19791 and dense parts of clouds become dimmer in the CO 


lines. 

Note that in our simulations gas number density can 
be high as 2000-3000 cm“®. However, even in such dense 
medium a star does not form with necessity, because a gas 
can be in the equilibrium with the surrounding medium. 
Such picture is usually taken place in small clouds. So that 
sometimes one can find rather small clouds (see Figs. |4f,b) 
with large amount of molecular gas, and these clouds ap¬ 
pear to be brighter than it is expected from third Larson’s 
relation. Thus, in our calculations the clouds surface den¬ 
sity is expected to be not always constant that reflects that 
in our model there is no gas density threshold for the star 
formation process. 


5.4 Variation of the spectra resolution 

Spatial resolution in numerical simulation plays a signifi¬ 
cant role in understanding of the internal properties and 


basic physical parameters of GMCs. Fujimoto et al. (20141 


reported that the variation of spatial resolution strongly af¬ 
fects on the properties of cloud populations. At the same 
time, results of the PPV data cube analysis can depend on 
spectral resolution. In our previous simulations of synthetic 
spectra the velocity resolution equals Sv = 0.5 km s“^ , 
which is quite high for extragalactic observations. Although 
this value is comparable to that used in several studies (e.g.. 
Tan et al.|2013 1, most of the recent extragalactic surveys in 
molecular lines have been done with much lower spectral res¬ 
olution (Engargiola et al.|2003 Donovan Meyer et al.|20T3 


Schinnerer et al.|2013 Combes et al.|2014|. 

To check whether spectral resolution affects on the scal¬ 
ing relations, we calculate and analyze the PPV data with 
lower spectral resolution Sv = 1 and 5 km s“^ for the model 
F. In Tab we show the power-law indices for the scaling 
relations with different resolution values, we also combine 
the indices obtained in several observations with brightness 


temperature threshold equal to 1 K. We argue that the no¬ 
ticeable variations of the indices with Sv are due to that 
for lower spectral resolution small clouds are combined into 
larger ones on the line of sight when their relative motion 
and velocity dispersion is lower or comparable with the spec¬ 
tral resolution. Here we only report that there is a depen¬ 
dence of the clond population characteristics on spectral res¬ 
olution. An accurate quantitative consideration of this effect 
requires further detailed study. 


5.5 Variation of the threshold value 

In the previous subsections we have established that the scal¬ 
ing relations for the cloud ensembles obtained in our simula¬ 
tions are quite similar to these found in the observations. It 
is interesting to study a dependence of the power-law indices 
of the relations on value of threshold. 

We consider the relations obtained by both methods de¬ 
scribed in Section 3. Here we constrain ourselves to analyse 
the model of galaxy with the prominent structure - model F. 
To do this we vary threshold values in the following ranges: 
Vtot = (0.5 - 4) X 10^2 cm-2 for CDN and = (1 - 15) K 
for CF. Using the lower limits we extract clouds with mass 
less than 10® M© , while for the upper values at least 100 
molecular clouds remain in the catalogue. 

Figure 1^ presents the dependence of the power-law in¬ 
dices for three scaling relations on total column density 
threshold Nt^t (left panels) and brightness temperature 
threshold (right panels). The error bars correspond to the 
data dispersion obtained in the procedure of the 

power-law indices for a given threshold. Note that the num¬ 
ber of clouds definitely depends on threshold value, but it 
remains above 100 clouds in order to provide enough objects 
for statistics. 

For low Ntot we extract both extremely large and small 
clouds. Large clouds consist of a group of small clouds en¬ 
closed by the extended diffuse structure, which can be called 
as the intercloud medium. Such a structure includes both 
molecular and atomic gas. The increase of the threshold ex¬ 
cludes the intercloud medium. So that for higher threshold 
values we extract bright cores of the virialized clouds. One 
can see in Fig. (left panels) that the indices /3i for the 
simulated clouds change dramatically in the case of using 
the CDN method. Better agreement with the observational 
data can be found only with relatively low threshold values 
« 0.5 - 1 X 10^^ cm"^ 

Whereas the use of the CF method does not provide 
any signihcant variations of the scaling relation indices with 
threshold (see right panels in Fig. [^. In this case the cloud 
populations are vanished from capture of the intercloud 
medium, because this extraction method directly relates to 
regions with high molecular fraction. This explains that the 
result remains rather robust relative to variation of bright¬ 
ness temperature threshold value. The indices obtained by 
using the CF criterion are close to the observed ones in wide 
range of threshold values. Moreover, one cannot see signif¬ 
icant dependence on velocity resolution for the hrst scaling 
relation. For the other relations the decrease of resolution 
leads to the systematic shift of the index values, so that the 
dependence on threshold value remains more or less flat. 
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Figure 8 . The power-law indices for the scaling relations, 0i, i = 
1,2,3, as function of the total column density threshold for 
clouds extracted using CDN approach (left panels). In the right 
panels the indices Pi are shown as a function of the brightness 
temperature TjJ*' for the cloud samples found in analysis of CO 
line spectra using the CF method for various spectral resolution 
5v. 





5.6 Clouds mass spectra 

The indices of the GMCs scaling relations slightly vary for 
galaxies with different morphology, although both methods 
of the cloud definition are suffered from so-called environ¬ 
mental effects (see Figs. [^ [^ [^. In our case such effects 
come from remarkable large scale structures like spiral arms 
and galactic bar. 

To check the impact of the galactic environment on the 
clouds properties we calculate the cumulative mass functions 
for three types of galaxies, i.e. number of clouds N with 
masses Md greater than a reference mass M'j: 

Ar(M'i) = Af(Md < M'l). (9) 


The mass spectrum of molecular clouds is usually ex¬ 
pressed as a power law function (see e.g. Rosolowsky||2007 


Gratier et al. 20121, however more accurate approach is 


based on the truncated power law shape (Williams & Me 
Kee|1997 Colombo et al.|2014 1 which can be written in the 
form: 


iV(Mci < M'l) = No 


(^) -1 


( 10 ) 


where the index 7 shows how the mass is distributed among 
the cloud population. 

We compute the fits of the cumulative mass distribu¬ 
tion in the form ( |10[ ) (Fig.[^. For all considered models the 
slope of the mass distribution function 7 (Eq. 10 1 is greater 
than —2 that means that large massive clouds dominate in 
the total GMCs mass budget. One can see that the clouds in 
the CDN samples have rather steeper mass distribution than 
that for the CF sample. This demonstrates that most molec¬ 
ular mass tends to be concentrated in less massive clouds in 
the CF sample than that in the CDN one. In other words, 
small clouds are more numerous in the CF sample. That is 
clearly seen from Fig. and even from Fig. if one men¬ 
tion that the total masses of extracted clouds obtained using 
both methods are very close to each other. Such conclusion 
is general for all three galaxy models (see Fig. [^. 

A remarkable truncation of the mass distribution is seen 
for the CDN cloud samples in all models. This can be ex¬ 
plained by the engaging of numerous structures above the 
column density threshold in the dense environment. This 
suggestion is conhrmed by Fig. where one can see that this 
effect is more clear in the galaxy with the prominent spiral 
pattern and bar (model F) and weaker in the galaxy with¬ 
out structure (model B). Such significant truncation does 
not reflects the physical state of isolated molecular clouds, 
because the truncation is not detected for the CF sample of 
GMCs. 


As it was mentioned above (see Sections 5.1 5.2 5.31, 


there are not strong variations of the scaling relation indices 
on galactic morphology for the same CD criterion. However, 
the impact of the galactic environment on the GMCs prop¬ 
erties reveals distinctly. It is seen in the mass distribution 
profiles in Fig.|^ The shapes of the distributions for the CF 
samples are quite similar to each other: the 7 values are 
in the range [—1.70; —1.13], that once again shows a homo¬ 
geneity of these cloud samples. Especially this is remarkable 
for the galaxies without large-scale pattern: models B and 
H. Note that for the CDN samples the distributions coin¬ 
cide for M Pj 10® Mq . The mass spectrum in the model 
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Figure 9. The mass spectrum of GMCs extracted by using the 
CDN (three top lines) and the CF (three bottom lines) methods 
in the models of galaxies (Tab. The index 7 for the mass 
spectrum in the form is shown in the legend . 


F systematically differs from the others. We suggest that 
stronger stellar feedback and compression of GMCs in spi¬ 
ral arms taken place in the model F substantially affects 
on the GMCs mass distribution. However, the conformity 
of the scaling relations (see Table denotes that GMCs 
save their internal structure or the CD methods work in 
resemling way and as a result extracted structures (clouds) 
have rather close physical parameters. Similar influence of 
large scale structures can be noticed both in numerical sim¬ 
ulations of M83 ( Fujimoto et al.||2014 l and in observations 
of M51 (Colombo et al. 20141. More detailed discussion of 
such influence on statistical properties of clouds from the 


observational point of view can be found in (Hughes et al. 


2013). 


6 CONCLUSIONS 

In this paper we have presented a set of the galactic scale 
simulations of the Milky Way size galaxies of different mor¬ 
phological type: a galaxy without prominent structure, a spi¬ 
ral barred galaxy and a galaxy with flocculent structure. In 
our models we have taken into account star formation, stellar 
feedback, UV radiation transfer and non-equilibrium chemi¬ 
cal kinetics for CO and H2 molecules. Here we have focused 
on the statistical properties of molecular clouds obtained 
by two different extraction methods of gaseous structures. 
The first uses the total hydrogen column density threshold 
as a marker of the cloud border. The other cloud defini¬ 
tion method is based on extraction from position-position- 
velocity (PPV) data cubes for (1-0) line. Using both 

methods we have studied the empirical scaling relations 
known as Larson’s laws: ’velocity dispersion - cloud size’, ’lu¬ 
minosity - cloud size’ and ’virial mass - luminosity’ relations. 


Using our simulations we have created the position-position- 
velocity data cubes for several values of velocity resolution 
and have investigated how the physical parameters of clouds 
and the indices of the scaling relations depend on spectral 
resolution. Our results can be summarized as follows: 


• the number of spatially resolved molecular clouds in 
the simulations slightly depends on galactic type and equals 
~ 10®; size, mass, luminosity and other physical properties 
of giant molecular clouds obtained in the simulations are 
close to the observed ones in our and nearby disk galaxies; 
note that the physical parameters of clouds depend on the 
cloud definition method (see Figs. |3|5|6|7| ). 

• the diffuse (intercloud) gas can be catched using to¬ 
tal column density as threshold in the extraction of clouds, 
especially this can be significant in dense large scale struc¬ 
tures, e.g. within spiral arms or bar; such diffuse gas has 
higher velocity dispersion and lower CO line brightness in 
comparison with other cloud material, so that using this 
method of extraction we can not exclude overestimation of 
the ID velocity dispersion due to the projection effects even 
at high total column density threshold values; 

• giant molecular clouds found by using the 
CLUMPFIND (CF for shortness) algorithm have smaller 
sizes, masses and velocity dispersion than those extracted 
by using total column density as threshold. However, 
the distributions of virial parameter for both extraction 
methods show similar behaviour (see Fig. [^; 

• numerical models of galaxies with different morphol¬ 
ogy produce a substantial number of rather small GMCs 
{Rci < 20 pc) which are detectable by various methods 
considered in the paper (see Fig. &)■ This is more clear 
for position-position-velocity analysis, where large clouds 
are splited into smaller ones due to complex kinematics of 
gaseous flows. However the analysis of the mass distribution 
functions shows that the mass of the cloudy phase in galax¬ 
ies simulated here is mostly concentrated in large massive 
clouds (see Fig. [^; 

• physical parameters of GMCs depend weakly on galac¬ 
tic structure, namely, mass, size, luminosity and velocity dis¬ 
persion are locked in the same ranges for models of galaxies 
without structure, with prominent spiral pattern and with 
flocculent patten (see Figs. |5|6|7| ); indeed, we do not see sta¬ 
tistically sensible variations of the scaling relations in the 
models of galaxies with different morphology for a given CD 
criterion (see Table |^; however so-called environmental ef¬ 
fects can be clearly seen in the distributions of cloud masses: 
the mass spectra are steeper in galaxy with the prominent 
structure (see Fig. [^. 


Thus, we conclude that it is impossible to extract equiv¬ 
alent clouds populations by using two different clouds ex¬ 
traction methods considered here: the first is based on to¬ 
tal column density as threshold and the second is utilized 
the PPV data analysis. Obviously, the comparison between 
observational and simulated properties of GMCs should be 
based on the same cloud extraction technique. A significant 
role of the cloud definition method and selection criteria (e.g. 
spectral resolution, threshold value etc.) can correspond to a 
fact that the observable scaling relations for external galax¬ 
ies might not completely reflect real physical parameters of 
the ISM cold phase. 
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